NPS-S8DW9071A 


UNITED STATES 
NAVAL POSTGRADUATE SCHOOL 


AD 696.56 7 


| ‘i " + | 


DIGITAL ANALYSIS OF TURBULENGE DATA 
ON THE IBM 360/67 AT TEE 
NAVAL POSTGRADUATE SCHOOL 


by 


J. R. Wilson 
N. E. J. Boston 
W. W. Denner 


1 July 1969 


blis 
document has been arvroved for pu ; 
aieies and pals; its aistribution 15 unlisiteds 


Repro juced by the : 
PCLEARINGHOUS® 
for Federal Scient.:1¢ % Technical 

Information Springfield Va 22151 


eign. 
eee 2 


Siig RENAN  anestee te Dineter Fe cee Fe 


NAVAL FOSTGRADUATE SCHOOL 
Monterey, Califorma 


Rear Admiral Robert W. McNitc, USN R. F. Rinehart 
Suverintendent Academic Dean 
ABSTRACT: 


A system of time series programs used by the Institute of Oceanography 
of the University of British Columbia was made available to the Departmer* 
of Oceanography ef the Naval Postgraduate School in February 1969, This 
report summarizes the system and outlines the procedures to be followed 
in using the programs, 


The system consists of three programs labelled UBC FTOR, UBC SCOR 
and UBC FCPLOT. The program UBC FTOX% computes Fourier coefficients 
from selected channels of analog-to-digital tape and writes them on ar.- 
other tape. The program UBC SCOR reads the tape produced by UBC FIOR 
and from the Fourier coefficients calculates spectra, cospectra and quad- 
rature spectra for the channels indicated. These are computed for each 
data block. The printed output gives for each quantity the average, stan- 
dard deviation and a number representing the trend over the blocks. I[n 
the case of co~ and quad-spectra phase and coherence are also printed 
out. The program UBC FCPLOT provides a Calcomp plot of the spectra 
for qualitative analysis. 


These programs have been tested on the IBM 360/67 of the Naval 
Postgraduate School and produced for a test tape the same answers as 
produced by the U.B.C. machine. 


A system to develop the capability to use the SDS-9300 and the asso- 
ciated analog computer available at the Naval Postgraduate School to 
digitize data to be analyzed by the time series programs is included as 
Appendix I, 
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FORWARD 

The digital programs discussed in this report were developed at the 
University of British Columbia by John Garrett and Ron Wilson during 
1967. They have been applied to a large variety of geophysical turbu- 
lence data and found extremely successful. Because of the nature of the 
air/sea interaction program being carried out by the Department of Ocean- 
ography of the Naval Postgraduate School, such programs were desirable 
to be available at the computing facilities of the School. 

In February of 1969, Mr. Ron Wilson was hired as a consultart to 
convert the University of British Columbia programs for use on the IBM 
360/67 computer of the Naval Postgraduate School. This report is a sum- 
mary of the system based on notes provided by Mr. Wilson. The principal 
investigators take full responsibility for any errors or omissions that may 


occur in this report. 


Noe! Boston 
W. W. Denner 
Principal Investigators 


ABSTRACT 
: A system of time series programs used by the Institute of Oceanoaraphy 
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of the University of British Columbia was made availaLie to the Depart-~ 
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ment of Oceanography of the Nava! Postgraduate School in February 1969. 
This report summarizes the system and outlines the procedures to be fol~ 
lowed in using the programs. 

The system consists of three programs labelled UBC FTOR, UBC SCOR 
and UBC FCPLOT. The program UBC FTOR computes Fourier coefficients 
from selected channels of analog-tc-digital tape and writes them on another 
tape. The program UBC SCOR reads the tape produced by UBC FTOR and 


from the Fourier coefficients calculates spectra, cocpectra and quadrature 


spectra for the channels indicated. These are computed fcr each data 
block. The printed output gives for each quantity the average, standard 
deviation and a number representing the trend over the blocks. In the 
case of co~ and quad-spectra . phase and coherence are also printed out. 
The program UBC FCPLOT provides a Calcomp plot of the spectra for quaii- 


tative analysis. 
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to These programs have been tested on the IBM 360/67 of the Naval 
| Fostgraduate School and produced for a tast tape the same answers as 
produced by the U. B. C. machine. 
A system to develop the capability to use the SDS-9300 and the asso- ‘ 
ciated analog computer available at the Naval Postgraduate School to 
digitize data to be analyzed by the time series programs is included as 
Appendix I. 
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I. GENERAL 

This report is intended (a) to summarize the system of time series 
programs used by the Institute of Oceanography of the University of British 
Columbia which are now available to the Department of Oceanography of the 
Naval Postgraduate School and (b) to describe briefly the procedures to be 
folicwed in using them. 

These programs were made available in February 1969 and this report 
describes the programs as they existed at that time. Some updating and 
expansion of these programs may occur from time to time. A check on their 
current state is recommended prior to using. 

The turbulence data used to test these programs included downstream 
velocity fluctuations measured with a hot-wire anemometer (DISA U) and 
vertical velocity fluctuations measured with a sonic anemometer (Sonic W). 
These measurements were made simultaneously during August 1966 at the 


U.B.C. Institute of Oceanography experimental site at Spanish Banks, 


British Columbia. 


II DIGITAL SAMPLING PROCEDURE 


The system is based on the assumption that the data to be analyzed 
are initially recorded on analog tape. These analog data must then be 
converted to digital form before they can be processed by a digital com~ 
puter. The analog to digital procedure is of central importance and must 
be done with extreme care. Almost alway 3 there must be some conditioning 
of the signal prior to the analog-to~digital (A/D) conversion, Common 
problems to consider are analog signal level, aliasing and sampling fre- 
quency. Of course each analysis problem presents its unique difficulties 
and must be treated individually. 

A. Analog Signal Level 

The patch of signal to be analyzed may not have been recorded under 
ideal conditions. It may vary from almost too small (approaching system 
noise level) to almost too large (occasional clipping). To insure acod 
digitizing some amplification may be necessary t provide the optimum 
signal input to the A/D converter. If there is an appreciable DC level 
on the signal, high pass filtering or rerecording may be necessary before 
digitization can begin. In brief, the signal level must always be checked 
to be sure matzimum advantage is being taken of the full dynamic range of 
the A/D converter. 

B. Aliasing 

When an analog signal is sampled at a frequency F, energy present 
in the signal at frequencies f >t will appear in the sampled output at an 
apparent (aliased) frequency F-f. This phenomenon ir known as aliasing 
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with called the folding or Nyquist frequency (Blackman & Tukey, 1951). 
To prevent aliasing, the input to the analog-to-~digital converter must 
contain no energy at frequencies greater than the folding frequency asso- 
ciated with the particular sampling rate being used. This can be done by 
i using an analog low pass filter between the signal source and the converter 
input. The filter cutoff frequency is normaily less than the folding fre- 
quency. In practice the’ is no exact relation between these frequencies 
as it depends on the input levels expected at high frequencies and the 
sharpness of the filter cutoff. 
C. Sampling Rate 
If three channels are to be sampled and co~spectra are required 
between all combinations the sampling procedure might be as iliustrated 


in Figure 1. 


Channel 1 


Channel 2 


Channel 3 


Cl con 
ccsl CCS2 
————~ ] cycle of highest frequency of interest ~————» 


Figure 1. Sampling Procedure 


10 


The cuts Cl, C2 and C3 represent the digital samplings of channels 
1, 2 and 3 respectively. 

CCS1 and CCS2 represent two samplings of all channels and are 
referred to as Cross Channel Sweeps. 


The time lag between sampling channel 1 and channel 2 in a crass 


: e 
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channel sweep should be a small fraction of a cycle since the program 
assumes they are sampled simultaneously for purposes of computing the 
co-spectra and quadrature spectra. 

The time between cross channel sweeps should be such that each 
channel is sampled at 1east twice per cycle of the highest frequency with 
sufficient amplitude to produce an aliasing problem. 

D. Miscellaneous Factors 
It is important that all factors influencing the validity and interpre- 


tations of results be considered prior to digitizing. Besides those factors 
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already discussed these should include the record length, the block length, 
spectral windows and how these will effect the stability of the spectral 
estimates. Further, an estimate of the computer time required for the 
analysis should be made. There has not been sufficient experience with 
the computing facilities at the Naval Postgraduate School to allow tables 

of time for typical calculations to be presented. However the tables given 


by Garrett (1967) will indicate order of magnitude computing times. In 


addition his report contains excellent discussions of some of the factors 
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mentioned here. 
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ill COMPUTATION OF SPECTRA AND CROSS SPECTRA 

The computation of the spectral densities employ the Fast Fourier E 
transform (FFT) of Cooley and Tukey (1965). The program used is dis- - 
tributed through ti. share library and is catalogued at PKFORT SDA3465. 

The subroutine will transform a number of data points equal to 9h 
where n is a positive integer less than or equal to 13. This means the 
largest number of data points that can be transformed is 8192. Since the 
time series encountered in turbulence measurements generate more points 
than this, some method of transforming and averaging consecutive blocks 
of data must be used. This also makes it possible to gain some information 
as to the variance of the spectral or cospectral densities over the interval 
of the signal. 

The method employed consists of evaluating the spectrum, co- 
spectrum or quadrature spectrum for each frequency band for successive 
blocks of 2” data points. The megan and variance of these densities are 
then formed and printed and plotted. 

For example purposes, assume 100 blocks of 2048 points each are 
transformed. For each biock there will be contained in the Fourier coef- 
ficients an amplitude and phase for frequencies corresponding to 
0,1, 2, ..., 1024 cycles over the 2048 points. Zero is the mean and 
1024 is the folding frequency. Therefore for each block of 2048 data points 
there are 1024 estimates of spectral, cospectral or quadrature spectral 
density. This is considered to be too much information to plot graphically 
so in the production of printouts and plots further averaging than along 
block must be done. For example the program may divide the spectrum 
into 32 bands and average across these bands as well as along blucks to 


obtain the mean and variance of the quantities. 
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If there are two channels and they have for the jth harmonic, ampli- 


tudes Ai and Bi , and phases a, and B, , then the spectral densitie + 
are Aa 32 


i 


and 2 x Bandwidth 


i 
2 x Bandwidth 
where the bandwidth is given by 


sampling frequency 
no. of data points in the trans‘orm 


The cospectrum is given by 


Ay x B. cos (a, - 3) 
2 x Bandwidth 


The quadrature spectrum is given by 
A, xB, sin (a, -3;) 
2 x Bandwidth 

These quantities are computed for each frequency and block anc are 
then averaged across the frequencies and along blocks. 

The error bars plotted are computed by evaluating the variar.ce ot the 
estimates of the spectra’ density in the band and dividing by the square 
root of the number of estimates to establish the expected variance in the 
mean. In the case of the SCOR program the plotted and printed error bars 
represent standard deviations in the mean. The FCPLOT program uses 
1.96 times the standard deviation (or some other kind of 95% confidence 
interval). 

The values of phase and coherence produced by the SCOR program 
are erived directly from the final values of the co- and quadrature spectra 
appearing to their left (Table II', for example). 


Further information is available on the comments cards with the 


Fortran source programs ‘Appendix II). 
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: Examples of print out for spectra, co-spectrum, quadrature spectrum 


and Fourier coefficients are given inTables I through V. The corresponding 


FCPLOTs for these tables are Figures 2 through 6. 
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IV INPUT DATA TAPE FORMAT 

The input data tape is read at the Fortran level by an unformated 
read. The read is in the subroutine OCEAN 1 and has the form: 

READ (NUNFT) ICTR, JCHANS, (DATA(J), J = 1, JCTR) 
where JCHANS is the number of channels digitized and JCTR is the number 
of data words in the tape record. DATA is an array into which the data is 
to be read and is single precision real and dimensioned 256. 

To create a suitable tape with 4 channels of information the programmer 
would set JCHANS = 4 then line un-floating point data words in a single 


dimensioned array, DATA, in the order: 


[ Channel 1 DATA(1) 
Cross Channel 2 DATA(2) 
Channel [ome 3 " (3) 
Sweep 1 Channel 4 " (4) 
Channel 1 " (5) 
Cross Channel 2 " (6) 
Channel Channel 2 " (7) 
Sweep 2 Channel 4 " (8) 
_ " (9) 
| Channel 1 
| Cross Channel 2 
| Channel Channel 3 ‘ 
| Sweep K Channel 4 DATA(JCTR) 
| 


' where 


ar JCTR=4*K and is less than 256. 
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The statement 


WRITE(KUNF7Z) JCTR, JCHANS, (DATA(J), J = 1, JCTR) 
ts then executed. 

This loading of the array and writing it out is repeated until all data 
has been transferred to tape and then an end of file is written. It is not 


necessarv thet JCTR be the same for all records written and it may be 


zero. It must be however an integral number of cross channel sweeps. 


V PROCEDURE IN ANALYZING DATA 
The first requirement is to produce a tape containing the data in the 
format that has been described. This tape is used as input to the FTOR 
program and a second tape is produced containing the Fourier coefficiterts. 5 
This tape is then used with either the SCOR or FCPLOT program depending on 
the analysis desired. If co-spectra or quadrature spectra are desired then 
SCOR must be used. If power spectra alone are desired SCOR or FCPLOT 
can be used. If more detail than SCOR supplies for the spectra, ora 
calcomp or printed plot of the individual Fourier coefficient amplitudes are 
desired, FCPLOT should be used. 
The format of the deck. of control cards required to run each of the 
programs is described in detail on the comments cards with the Fortran ie 
source decks (APPENDIX) 
The programs were originally written for an IBM 7044 (Garrett, 1967) 
and were then modified for an IBM 360/67 running under the Michigan 
Terminal System. Neither of these systems have the set of utilities for the 
tape handiing available under OS/MUT. For this reason all tape positioning 
on files is handled by the programs at the Fortran ievel. It is therefore 


necessary to have a DD card present in the deck for each file on the tape 


whether it is read or written or just spaced over, 


SES eae tego mn 2 ARR OSM RTS OST RTPI. MEP MN SE SRN ET aa 


VI SOURCE PROGRAMS 
The system consists of three programs: 
1. UBC FTOR 
2. UBC SCOR - 150K in FORT step to compile. 
3. UBC FCPLOT 
The source programs have been written as files 1, 2 and 3 in the order 
above on the tape NPS216. The DSNAME's are UBCFTOR, UBCSCOR sand 
UBCFCPLT. The DCB is (RECFM = FB, LRECL = 80, BLKSIZE = 1600). 


The object programs ere on disc pack FAC001 in load module form. The 


DSNAME for the file is F1178.TSLIB. 
UBCFTOR 
UBCSCOR 


UBCFCPLT 
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APPENDKX 1 
ANALYZING DATA DIGITIZED ON THE SDS-9300 


Two computer programs are required to develop the capability to use 
the SD$-9300 and the associated analog computer to digitize data to be 
analyzed by the time series programs. 

The first is an SDS-9300 program to control the digitization procedure 
and produce a 7-track magnetic tape containing the digital data. The 
second program is an IBM 360 program to convert the 7-track tape toa 
9~track tape in the format described in this report for the input tape to the 
time series program. 

The SDS-9300 program should do approximately the following: | 

1. 1tshould allow for manual starting and stupping of the digitization ". ! 

process by the computer operator when provided with suitable sig- | 


nals (audio tones, oscilloscope signal changes, verbal instructions 


from the person concerned with the analog signal, etc.) 
2. It should allow for a variable number of analog input channels. 

This number should be an input to the program at execution time. 
‘ Ten would be a safe maximum since the time series system of 
programs can handle only ten. 
| 3. The program should be written so that when the interrupt is given 
| all channels are sampled as rapidly as possible and the data accu- 
mulated for tape output. Only after a cross channel sequence is 
completed should a test and possible tape output occur if enough 


1 
data has been accumuiated. 


4. Care should be taken to make the program efficient in the sense 


that the computer is always ready to service the clock interrupts 
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when they cccur. Further, care should be taken to make the 
records long enough that the tape is used efficiently. 

The program for the IBM 360/57 should take this 7-track tape 
and produce a 9-track tape in the format required by the time 
series programs. 

It is anticipated that multifile tapes wiil be necessary to a'7oid 
piling up large numbers of magnetic tapes with only the first 
part of the tape used. Therefore these two programs should be 
made capable of positioning themselves in multifile tapes under 
the control of either the job control language or input control 


cards to the program. 
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APPENDIX I] 
FORTRAN SOURCE PROGRAMS 
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OCEAN3 
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This project funded Ly: Naval Postgraduate Schoo} 
Naval Ordnance Systems Command Monterey, California 


* ***"R*System of time series programs used by the Institute of Oceanography of the Uni- 
versity of British Columbia was made availabe to the Department of Oceanography of 
the Naval Postgraduate Schocl in February 1969. This report summarizes the system 2 
outiines the procedures to be followed in using the programs. 


The system consists of three programs labelled UBC FTOR, UBC SCOR and UBC FCPL@T. 
The program UBC FTOR computes Fourier coefficients from selected channels of analog- 7 
to-digital tape and writes them on another tape. The program UBC SCOR reads the tape 
‘ produced by UBC FIOR ani from the Fourier coefficients calculates spectra, cospectra 
and quadrature spectra for the channels indicated. These are computed for each data 
block. The printed output gives for each quantity the average, standard acviation and 
a number representing the trend over the blocks. In the case of co~ and quad-spectra 
phase and coherence are also printed out. The program UBC FCPLOT provides a Cal- 
comp plot of the spectra for qualitative analysis. ,, 


These programs have been tested on the IBM 250/67 of the Naval Postgraduate Schoo 
and produced for a test tape the same answers as produced by the U.B.C. machine. 


l 

| A system to develcp the capability to use the SDS-9300 and the associated analog 
pe computer available at the Naval Postgraduate School to digitize data to be analyzed by 
the time series programs is included as Appendix I. 
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